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ABSTRACT 

Supersonic random motions are observed in dark clouds and are traditionally interpreted as 
Alfven waves, but the possibility that these motions are super-Alfvenic has not been ruled out. 

In this work we report the results of numerical experiments in two opposite regimes; f3 ss 1 
and (3 <!C 1, where (3 is the ratio of gas pressure and magnetic pressure: (3 = P g /P m . 

Our results, combined with observational tests, show that the model with (3 ~ 1 is consistent 
with the observed properties of molecular clouds, while the model with [3 <C 1 has several 
properties that are in conflict with the observations. 

We also find that both the density and the magnetic field in molecular clouds may be very 
intermittent. The statistical distributions of magnetic field and gas density are related by a 
power law, with an index that decreases with time. Magnetically dominated cores form early in 
the evolution, while later on the intermittency in the density field wins out. 


Subject headings: turbulence - ISM: kinematics and dynamics- magnetic field 


1. Introduction 

The observation of supersonic motions in molecular clouds (eg Zuckerman & Palmer 1974) raised the 
question of how these motions could be supported (Norman & Silk 1980; Fleck 1981; Scalo & Pumphrey 
1982). Supersonic motions are expected to quickly dissipate their energy in highly radiative shocks, because 
of the very short cooling time of molecular gas or metal rich atomic gas (Mestel & Spitzer 1956; Spitzer 
1968; Goldreich & Kwan 1974). 

Strictly related was the issue of the support of molecular clouds (MCs) against gravitational collapse, 
since it was soon realized that the observed motions could not be understood as a gravitational collapse 
(Zuckerman & Evans 1974; Morris et al. 1974), although MCs contain many Jeans’ masses. 

Theoreticians therefore formulated the hypothesis that MCs were primarily magnetically supported 
(Mestel 1965, Strittmatter 1966; Parker 1973; Mouschovias 1976a,b; McKee & Zweibel 1995) and interpreted 
the observed motions as long-wavelength hydro-magnetic waves (Arons & Max 1975; Zweibel & Josafatsson 
1983; Elmegreen 1985, Falgarone & Puget 1986). It was also shown that the properties of the observed 
‘turbulence’ (Larson 1981; Leung, Kutner & Mead 1982; Myers 1983; Quiroga 1983; Sanders, Scoville 
& Solomon 1985; Goldsmith & Arquilla 1985; Dame et al. 1986; Falgarone & Perault 1987) could be 
understood if the motions were sub-Alfvenic (Myers & Goodman 1988; Mouschovias & Psaltis 1995, Xie 
1996). 
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Nevertheless, recent attempts to detect the Zeeman effect, in lines of molecules such as OH and CN, 
that probe regions of dense gas (Crutcher et al. 1993; Crutcher et al. 1996), resulted in a number of 
non-detections of the effect, and therefore in rather stringent upper limits for the magnetic field strength, 
despite the fact that regions expected to favor detections were targeted. It is therefore possible, that the 
field strength in dense molecular gas is weaker than assumed in theoretical studies. 

In the present work, we describe the dynamics of MCs, with the numerical solution of the equations of 
3-D compressible magneto-hydrodynamics (MHD), in a regime of highly supersonic random motions. We 
avoid any sort of idealization of the physical system, apart from describing it as a fluid, and simply rely 
on the physical assumption (supported by the observations) that the motions in MCs are random. The 
main limitation of our numerical models, compared with MCs, is the absence of gravity. We have excluded 
gravity on purpose (our code is capable of handling self-gravity), because one of the aim of our work is to 
show that (magneto-)hyclrodynamic processes alone are able to explain most of the observed properties of 
MCs, without the need of gravity. Although gravity is certainly responsible for the fragmentation into stars 
of high density regions, supersonic random motions shape the structure of molecular clouds with a minor 
contribution from gravity (Padoan 1995; Padoan, Jones & Nordlund 1997a; Padoan, Nordlund & Jones 
1997b), apart from the possibility that gravity is the energy source of the motions on the large scale. The 
importance of supersonic motions in fragmenting the gas is certainly apparent on very small scales, where 
young and probably transient clumps are not originated by gravitational instability (Falgarone, Puget & 
Perault 1992; Langer et al. 1995). 

Previous numerical studies have shown that compressible turbulence can qualitatively explain several 
observational properties of MCs, even if the effect of gravity and magnetic fields are not included. The first 
two-dimensional (2-D) simulations of turbulence, with rms Mach number larger than one, were performed 
by Passot & Pouquet (1987). These were the first simulations where shocks are shown to develop inside 
a turbulent flow. The importance of shock formation inside mildly supersonic flows was later confirmed 
in 3-D simulations by Lee, Lele & Moin (1991). It was immediately recognized that shocks might have 
been responsible for the fragmentation of the density field inside MCs, and especially for the origin of their 
filamentary structure (Passot, Pouquet & Woodward 1988). Kimura & Tosa (1993) simulated the passage of 
a strong shock through a turbulent molecular cloud, and found that this process can generate dense clumps, 
with a power law mass spectrum. Vazquez-Semadeni (1994) made use of 2-D numerical simulations to show 
that supersonic turbulence generates a very intermittent density field, reminiscent of the clumpy nature 
of MCs. The density field was also found to be self-similar, which could be the reason for the hierarchical 
structure of MCs (Scalo 1985, Houlahan & Scalo 1992). Falgarone et al. (1994), analyzing the numerical 
simulation by Porter, Pouquet & Woodward (1994), argued that the properties of the profiles of molecular 
emission spectra from MCs can be interpreted as arising from turbulent motions. 

Other numerical works included gravity in the turbulent flows, yet without describing the magnetic 
field. Turbulence was shown to be able to prevent gravitational collapse (cf Chandrasekhar 1958; Arny 
1971; Bonazzola et al. 1987, 1992) in the 2-D numerical simulations by Leorat, Passot & Pouquet (1990). 
Vazquez-Semadeni, Passot & Pouquet (1995) modeled the galactic disc, on the scale of 1 Kpc, as a 
turbulent self-gravitating flow. They simulate a 2-D turbulent flow that is forced by the energy released 
by star formation (expansions of HII regions), and found that the main mechanism of cloud formation is 
the turbulent ram pressure, rather than gravity. They were not able to form self gravitating clouds, due to 
limitations in the thermal modeling, and the consequent low density contrast. 


A 3-D description of a magnetized self-gravitating cloud was given by Carlberg & Pudritz (1990), 
and was used to simulate molecular emission spectra by Stenholm & Pudritz (1993). Carlberg & Pudritz 
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found that the magnetic field and hydro-magnetic waves can support the cloud against gravity. The 
clouds contract, because of ambipolar diffusion, on a timescale of approximately four free-fall times. These 
simulations do not solve the MHD equations, but instead make use of a ’sticky particles’ code. Energy is 
injected in the form of a spectrum of Alfven waves, and the outcome of the computation is dependent on 
the spectral index, that is a free parameter. This way of forcing the particles is rather unphysical, because 
an arbitrary spectrum of Alfven waves is imposed, instead of being obtained as a result of the simulated 
magneto-hydrodynamics. Passot, Vazquez-Semadeni & Pouquet (1995) introduced the magnetic field in 
their previous 2-D model for the galactic disc (Vazquez-Semadeni, Passot & Pouquet 1995), and obtained a 
flow with rough equipartition of kinetic and magnetic energy, probably in rough equipartition also with the 
mean thermal energy. The same simulation, and others with larger density contrast and resolution, have 
been studied by Vazquez-Semadeni, Ballesteros-Paredes & Rodriguez (1997), who were able to reproduce 
the observed relation between line-width and size (Larson 1981). Gammie and Ostriker (1996) solved the 
MHD equations in a slab geometry, including self-gravity. By forcing the flow with a nonlinear spectrum of 
MHD waves, they were able to prevent the gravitational collapse of their 1-D cloud model. 

Apart from the intentional exclusion of gravity, we have improved significantly on all previous 
calculations of turbulent flows. First of all we have solved the MHD equations in three-dimensions, while all 
previous solutions are in two or one dimensions. In MHD the dimensionality of the flow has a fundamental 
importance, because it determines the topological freedom of the magnetic field. Another improvement of 
our work is the high rms Mach number of the flows (~ 5), while previous models are only mildly supersonic 
(Mach number ~ 1). While previous models of magnetized clouds focused on the role of the magnetic field 
as opposed to gravity, and therefore have assumed a magnetic pressure much larger than the gas pressure, 
the main purpose of the present work is to show that MCs are well described as flows with much lower 
magnetic pressure than previously assumed, and probably in rough equipartition with the gas pressure. 

We report on the results of two numerical models. In model A, (3 « 1, in model B, (3 <C 1, where (3 is 
the ratio of gas pressure and magnetic pressure: (3 = P g /Pm- The observations of magnetic field strengths 
are in agreement with a scenario where the mean magnetic pressure, ( B 2 a ), is dynamically low (model 
A). Magneto-hydrodynamic (MHD) flows develop a very intermittent spatial distribution of the magnetic 
energy, and therefore, when the field is detected at a favorable position, its strength is far stronger than the 
mean field strength. Dense cores with sub-Alfvenic velocity dispersion can still be generated, in agreement 
with the observations. 


2. The experiments 

The study of the dynamics of MCs belongs to the field of random MHD flows, in a highly supersonic 
regime. The Reynolds number and the magnetic Reynolds number in MC flows are very large. Their 
random nature is therefore a basic feature of the dynamics of these flows, and requires an appropriate 
description. 

For this reason, a realistic description of the dynamics of molecular clouds had to be based on the 
numerical solution of the compressible MHD equations in three dimensions, in a random and highly 
supersonic regime. We did it with relatively high resolution (128 3 ), with a code designed for turbulence 
and MHD turbulence experiments (Nordlund, Galsgaard & Stein 1994; Stein, Galsgaard & Nordlund 1994; 
Galsgaard & Nordlund 1996; Nordlund, Stein & Galsgaard 1996; Nordlund & Galsgaard 1997), specifically 
adapted to be able to deal with very strong shocks and very large density contrasts. 
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One of our purposes with these particular experiments is to explicitly address the classical arguments 
about dissipation time scales for super- and sub-Alfvenic motions, and we therefore perform the experiments 
in the spirit of “freely decaying turbulence”; studying the decay of an initial, solenoiclal velocity held, 
without applying any external forcing. We include external forcing in other experiments that we have made 
(e.g., Padoan et al. 1997a, 1997b). 

Although we have already developed a version of the code with the inclusion of the gravitational force, 
all the experiments were run without gravity, for the following reasons: 


• We are mainly interested in studying the magneto-hydrodynamics, rather than the gravitational 
instability. 

• The observed motions have velocities comparable with the virial velocity, or larger, on a range of 
scales, and the clouds are not free-falling. 

• If the results of our experiments are discussed only up to a time shorter than or comparable to the 
dynamical (or free-fall) time, all our conclusions remain basically unchanged. This time is about a few 
million years on a scale of 10 pc, and clouds are not supposed to live much longer than that, before 
star formation takes place and becomes energetically important. 


2.1. The equations 


We solve the compressible MHD equations: 
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( 4 ) 


j = VxB. (5) 

plus numerical diffusion terms, and with periodic boundary conditions, v is the velocity, B the magnetic 
field, f an external force (= 0 in these particular experiments), and p = pT is the pressure at T k const. 

Conditions in the cold molecular clouds that we are modeling are such that an isothermal approximation 
is adequate; the radiative heat exchange is so efficient that the temperature remains low in most places. 
Even if the temperature momentarily increases in shocks, the subsequent cooling is rapid, and the result is 
shock structures that are qualitatively and quantitatively similar to isothermal shocks. 
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We have thus used isothermal conditions in most of our runs and have verified that this is appropriate, 
by rerunning segments of some experiments using the full energy equation. No significant change of 
the statistics was found and, since using the full energy equation increases the cost of the experiments 
considerably (the strong cooling required to maintain a low temperature forces a much smaller time step), 
we performed most of the experiments at constant temperature. 

The absence of an explicit resistivity p in the induction equation corresponds to an assumption of flux 
freezing on well resolved scales. The code uses shock and current sheet capturing techniques to ensure 
that magnetic and viscous dissipation at the smallest resolved scales provide the necessary dissipation 
paths for magnetic and kinetic energy. As shown by Galsgaarcl & Nordlund (1976, 1977), dissipation of 
magnetic energy in highly turbulent MHD plasmas occurs at a rate that is independent of the details of 
the small scale dissipation. In ordinary hydrodynamic turbulence the corresponding property is one of the 
cornerstones of Kolmogorov (1941) scaling. 

We have not included ambipolar diffusion in the present experiments, because it occurs on a time-scale 
significantly longer than the dynamical time, as recently shown by Myers & Khersonsky (1995) and is thus 
expected to be of secondary importance. However, our code has the ability to handle ambipolar diffusion, 
and we plan to study its effects in a subsequent paper. 


2.2. The code 


The code solves the compressible MHD equations on a 3D staggered mesh, with volume centered 
mass density and thermal energy, face centered velocity and magnetic field components, and edge centered 
electric currents and electric fields (Nordlund, Stein & Galsgaard 1996). 

The original code works with “per-unit-volume” variables; mass density, momenta, and thermal energy 
per unit volume. In the super-sonic regime relevant in the present application, we found it advantageous to 
rewrite the code in terms of “per-unit-mass” variables; Inp, u, and E = With these variables, the time 
evolution of all variables is governed by equations of the type 


££ 

Dt 


9l 

dt 


^ J , £ 

— + v-/ = ...; 


( 6 ) 


i.e., equations that specify the time rate of change following the motion. These are better conditioned 
than the divergence type equations that result from using per-unit-volume variables (the large—order 
M 2 — density variation in isothermal shocks cause the per-unit-volume fluxes to vary over several orders of 
magnitude). 

We use spatial derivatives accurate to 6th order, interpolation accurate to 5th order, and Hyman’s 3rd 
order time stepping method (Hyman 1979). 

In order to minimize the viscous and resistive influence on well resolved scales, we use monotonic 
3rd order hyper-diffusive fluxes instead of normal diffusive fluxes, and in order to capture hydrodynamic 
and magneto-hydrodynamic shocks we add diffusivities proportional to the negative part of the velocity 
divergence, and resistivity proportional to the negative part of the cross-field (two-dimensional) velocity 
divergence. Further details of the numerical methods are given by Nordlund, Galsgaard & Stein (1996) and 
Nordlund & Galsgaard (1997). 
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2.3. Weak and strong magnetic field 

For the purpose of this work we have run two experiments: one with /?$ = 4 (model A), and the other 
with /3i = 0.04 (model B), where /% = (P g /P m )i is the initial ratio of gas and magnetic pressure. 

In both experiments the initial density is uniform, and the initial velocity is random. We generate 
the velocity held in Fourier space, and we give power, with a normal distribution, only to the Fourier 
components in the shell of wave-numbers 1 < k < 2. We perform a Helmholtz decomposition, and use 
only the solenoidal component of the initial velocity. The initial magnetic held is uniform, and is oriented 
parallel to the z axis: B = BqZ. 

In experiment B, the initial velocity held is perpendicular to the magnetic held (zero 2 component), in 
order to excite the Alfven waves, while in experiment A the initial velocity held has equal amplitude in all 
three components. Both experiments are decaying hows, because no external forcing is applied. 

We dehne the Alfvenic Mach number, A4 a, as the ratio of the how rms velocity and the Alfven velocity. 
In model A, A4a ~ 10, and in model B, A4a ~ 1, initially. Under isothermal conditions, in our units, the 
ordinary Mach number M. = A 4a/ P 1 ^ 2 , so the ordinary Mach number is initially A4 s=s 5 in both cases. 

Fig. 1 shows the time evolution of the rms Mach number and rms density in the two experiments. In 
both cases, the rms density grows to the value o n ss 0.5A4, where M is the rms Mach number of the how, 
in about one dynamical time, dehned as the ratio of half the linear size of the box over the initial rms 
velocity, t dyn = 0.5L box /a Vt0 - 

The initial Mach number in experiment B is a bit smaller than in experiment A, because the amplitude 
of the initial 2 component of the velocity held is zero. Nevertheless, energy is immediately transferred to 
the 2 component of the velocity, and the rms 2 component is subsequently about half the value of the other 
two components. 

Fig. 2 illustrates the time evolution of the kinetic and magnetic energies, expressed in units of the 
mean thermal energy. The time evolution of the magnetic energy is totally different in the two experiments. 
In model B the held lines are only weakly perturbed by the how, and exchange their energy with the how 
periodically (oscillations in magnetic and kinetic energy, E m and Ek)', in model A, instead, the magnetic 
held is almost passively advected, and its intensity and direction are strongly influenced by compressions 
and stretching of held lines. E m grows until it is in approximate equipartition with the thermal energy of 
the gas, but remains well below equipartition with the kinetic energy. 


3. Observational tests 

In this section we compare our numerical results with observed properties of molecular clouds. 


3.1. Dissipation of supersonic motions 

The time evolution of the kinetic energy is plotted in Fig. 2. After two dynamical times, t = 2.0 tdyn, 
the how in experiment A still contains about 30% of its initial kinetic energy or, in other words, the rms 
velocity is still about 55% of its initial value. In experiment B, the mean kinetic energy of the oscillations 
after two dynamical times is about 40% of its initial value, i.e. the rms velocity is about 63% of the initial 
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value. 

Thus, contrary to the beliefs that partly motivated developing the Alfven wave model of MC turbulence, 
the dissipation in model A is not particularly rapid in absolute terms, and the dissipation in model B is 
not particularly small. Thus, the advantage of the Alfven wave model, as far as the energy dissipation 
time-scale is concerned, is little: tdi S s,B = 1.3 tdi S s,A- Although this was already appreciated in Zweibel & 
Josafatsson (1983) and Elmegreen (1985), it should be noticed that here we do not include the dissipation 
mechanism of ion-neutral friction. Alfven waves dissipate mainly because they transfer their energy to 
motions along the magnetic held, that are eventually dissipated in shocks. This was found also by Gammie 
and Ostriker (1996). It is clear therefore that, from the point of view of the dissipation mechanism, model 
B is not significantly superior to model A. 

The dissipation of the supersonic motions is not a problem for any of the models, however, because: 

• If motions are present there must be a source for their energy; such a source may be active for more 
than one dynamical time. 

• After one dynamical time 50% of the initial energy in model A is still in the flow. At the scale of 10 
pc the dynamical time is about 3 x 10 6 years: clouds do not need to exist much longer than that 
before star formation significantly affects their energetics (Blitz & Shu 1980). 


3.2. Magnetic energy probability distribution 

We illustrate the probability distribution of the magnetic energy in units of the mean thermal pressure: 
1//3 = B 2 / < P g ). We recall that initially (3 is uniform in both models, and its value is = 4 in model A, 
and /3i = 0.04 in model B. 

The distribution of the magnetic energy is shown in Fig. 3. Model A develops a very intermittent 
distribution, with an exponential tail. 

Taking into account the relation between field strength and density (see next section), one finds that 
approximately 0.5% of the total mass of the system contains a field 10 times stronger than the mean value. 
For example, in a cloud of 10 4 Mq, one can find a couple of clumps of about 30Mq each, with a held 
strength B = 40/rG, while the mean held of the cloud is only B = 4fiG. 

Model B is much less intermittent than model A. 


3.3. The relation B-n 

In regions of maser emission, at densities of about n = 10'cm -3 , a held strength of the order of 
B = 10 3 — 10 4 ^G is observed, while in regions of molecular emission, with approximately n = 10 2 — 10 3 cm -3 , 
the held is found to be of the order of B = 10/iG (eg Myers & Goodman 1988). 

A relation of the type B oc n 0 3 ^ 0 6 may be deduced from the observations (eg Troland, Crutcher & 
Kazes 1986; Heiles 1987; Dudorov 1991), but it is quite uncertain, especially in the light of the above 
discussion about the intermittency of the distribution of the magnetic energy. 

The fact remains, however, that observationally the held strength certainly grows with the gas density, 



and this is found over a range of 6 orders of magnitudes in density, and 3 orders of magnitude in field 
strength. 

Fig. 4 shows the relation B — n in models A and B. The dispersion in the relation B — n is shown in 
the form of 1 — a ‘error’ bars in the plot of Fig. 4. In model A the slope of the relation decreases with 
time. We find B oc n 0 ' 8 at t = 0.2 td yn , B oc n 0 7 at t = 1.0 tdyn, B oc n 0 ' 5 at t = 2.0 tdyn, and B oc n 0 ' 3 at 
t = 3.0 tdyn- We interpret this time evolution of the relation B — n in the next section. 

In model B there is no significant correlation between B and n. The lack of correlation between B and 
n was to be expected: most of the density enhancement occurs by convergent flows along the field lines, in 
the 2 direction, and these motions do not affect the magnetic field. A nice illustration of this is provided 
by the snapshots in Fig. 5 (lower row), where the density field is stratified in planes that are predominantly 
perpendicular to the z direction (vertical direction in the images). The densest regions develop at nodes in 
the large-scale waves, as found by Carlberg & Pudritz (1990). 

One could possibly argue that the lack of correlation of B and n is consistent with the majority of 
observations, because no field is detected, and hence nothing can be said about the correlation. However, as 
further discussed in Section 4 below, the upper limits of the non-detections speak against this interpretation. 
Also, the cases where a field is detected would then have to be explained with ad hoc arguments, rather 
than as a natural part of a statistical distribution. 

An interesting consequence of the model A results is that the observed relation could in principle be 
used to set an independent estimate of the age of molecular clouds: the comparison between model A 
and the observations would indicate that molecular clouds, on the average, have survived, since they were 
almost uniform in density and field distribution, for 2 or 3 dynamical times, that is about 10' years, on a 
scale of 20 pc. This estimate has little to do with the estimate of the time-scale for massive stars to destroy 
the parent molecular cloud (Blitz & Shu 1980), since that calculation has nothing to say about the lifetime 
of the cloud after it cooled from diffuse warm gas. 

To conclude, note that Fig. 4 illustrates that, even if in model A the mean magnetic field is such that 
the (f3) « 0.7, in regions only 10 times denser than the mean, the field can be as intense as in model B. 


3.4. Cloud and flow topology 

An understanding of the spatial structure of the density field in dark clouds is very important for a 
correct interpretation of observational data, and for the formulation of a number of physical models. 

The snapshots in Fig. 5 give an idea of the dimensionality of the structures in the density field. It 
is apparent that experiment B (lower row of panels) develops a stratified density field, with sheet-like 
structures, the sheets being approximately perpendicular to the magnetic field (the z direction corresponds 
to the vertical direction in the figures). This is consistent with the fact that only the motions in the z 
direction can compress the gas to very high density. 

The topology of the density field in experiment A (upper row of panels in Fig. 5) has a clear evolution 
in time. In the very beginning, until t ss 0 Atd yn , the density grows predominantly in sheets. These are the 
fronts of blobs of coherent motion, advancing at supersonic velocity. Later, these fronts start to intersect 
each other, and the density increases especially in filaments (at the intersections of fronts). The evolution 
continues with the intersection of filaments into knots, at t ss 1.5 tdyn- The fully developed topology, at 
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t = 2.0td yn , is characterized by both filaments and knots (cores). 


The evolution of the magnitudes of the mass density and the magnetic field may be discussed with 
reference to Lagrangian version of the continuity equation, 


Dlnp 

Df 


-V-v, 


(7) 


and the scalar induction equation 


D In \B 
Dt 


— V_L ' V, 


( 8 ) 


where Vj_- stands for the divergence perpendicular to the magnetic field, following both the motion and the 
change of orientation of the held lines. 


Although — V • v vanishes for the solenoidal initial condition, the supersonic motions rapidly lead to the 
formation of shock fronts, where the local value of — V • v is large and positive because of the discontinuity 
in the velocity perpendicular to the shock. 


The initially homogeneous magnetic held is carried along by the perpendicular components of the 
velocity held, and is hence also collected into sheets, except at those rare locations where the initial held 
happens to be strictly parallel to the velocity held. This explains why sheets initially form in both mass 
density and magnetic hux density, and why the B — n relation initially has an exponent close to unity. 

Note that the usual argument for a slope of 2/3, that applies to isotropic and non-shocking compressible 
motion does not apply here, because of the development of discontinuities. In term of Eqs. 0 and ||, the 2/3 
follows if — V_l • v typically picks up two of the three directional contributions to —V • v. However, at a 
shock, the divergence is dominated by the derivative in one particular direction; the one perpendicular to 
the shock front. The magnetic held that is swept into the discontinuity quickly becomes almost parallel to 
the shock front, because the component in the plane of the shock grows exponentially with time. Thus, as 
long as the topology is dominated by sheets, the mass density and the magnetic hux grow more or less in 
unison in the sheets, corresponding to an exponent in the B — n relation close to unity. 


In the subsequent evolution, there are effects that tend to reduce the exponent in the B — n relation. 
First, the non-linear evolution of the initially solenoidal velocity held also leads to the development of 
regions of space with a positive divergence, in which both the mass density and the magnetic hux density 
decline. In these regions, there is no particular dominance of the cross-held divergence, and thus the 
three-dimensional divergence picks up un additional contribution relative to the two-dimensional divergence, 
consistent with the classical 2/3 argument outlined above. 

In experiment B the motions are mainly Alfven waves, and therefore the velocity is predominantly 
perpendicular to the direction of the magnetic held. This is illustrated in Fig. 6, where we have plotted the 
histogram of cos(a), where a is the angle between v and B. Note that, even though motions across the 
held lines are the most common ones, it is the less frequent motions along the held lines that dominate the 
dissipation, as mentioned in connection with the discussion of Fig. 5 above. The motions across held lines 
are subject to magnetic restoring forces, and do not lead to density enhancements. It is the motions along 
the held lines that lead to the sheet like density enhancements visible in Fig. 5. 

In experiment A, the magnetic held is advected by the how, and the stretching of held lines instead 
produces some alignment between v and B, already before one dynamical time has passed, as illustrated in 
Fig. 6. 

Alignment between B and v may be caused by two, complementary effects: 1) Dynamical alignment 
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is expected when the magnetic energy approaches and exceeds the kinetic energy; the Lorentz force then 
forces the flow to be predominantly along the magnetic held lines. 2) Kinematic alignment occurs when a 
spatially non-uniform velocity held causes stretching of magnetic held lines, and hence a correlation of B 
and v. Pure shear, for example, tends to align an embedded magnetic held with the direction of the how. 

Motions that are aligned with the magnetic held affect the mass density without affecting the magnetic 
hux density. In particular, the non-linear concentration into first sheets and then hlaments due to the 
interaction of shock fronts continues into the formation of knots in the density held, by the convergence of 
matter howing along hlaments. There is no corresponding process available to a divergence free vector held 
such as the magnetic held; once the held has concentrated into hlaments, it cannot concentrate further; the 
magnetic held in a filament is insensitive to how along the filament. 

In the same way that converging hows along the magnetic held may lead to extreme concentrations of 
mass, those regions where the how is diverging along magnetic held lines may lead to extreme rarefactions 
of mass, without affecting the magnetic hux density. In B — n scatter plots, this corresponds to the 
development of more extreme excursions of the mass, relative to those of the magnetic hux density, and 
hence a flattening of the B — n relation with time. 

In model A, dynamical alignment is at most significant in the few cores that develop a strong 
(sub-Alfvenic) magnetic held in the early evolution of the experiment. In scatter plots of B against n 
most contributions come from regions where dynamical alignment is unimportant. We thus conclude that 
the evolution of the B — n relation in model A, towards a smaller exponent with time, is caused by the 
kinematic alignment of B and v. 


3.5. Stellar extinction 

Padoan, Jones, & Nordluncl (1997) have shown that near-infrared stellar extinction determinations 
can be used to infer the 3-D probability distribution of the gas density in dark clouds. They have shown 
that there is qualitative and quantitative agreement between the inferred 3-D density distribution in dark 
clouds, and the one produced by their experiment of random supersonic hows. 

The method is based on the plot of the dispersion of the extinction measurements in cells, versus the 
mean extinction in the same cells (Lada et al. 1994). In Padoan, Jones & Nordlund (1997), the theoretical 
plots are produced by generating random density fields of given statistics and power spectra. Here we show 
one example of the same plot, but produced directly from the density field of experiments A and B. The 
plot is shown in Fig. 7, that can be compared directly with the plots shown in Padoan, Jones, & Nordlund. 
The plot of model A looks like the observational one. It has a smaller slope because the rms Mach number 
in Experiment A is smaller than in the observed cloud. 

The plot for Experiment B seems instead to differ from the observational one. The absence of locations 
with very low extinction in model B is due to the fact that the density field in experiment B is mainly 
organized in sheets, perpendicular to the magnetic field direction. 
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4. Discussion 

It is difficult to get an objective view on the magnetic field strength in dark clouds from the literature. 
The reasons are the following: 

• Negative results from observational programs, in which detections have been reported, often remain 
unpublished. 

• The positions searched for Zeeman splitting never represent a statistically meaningful sample. 
Favorable regions are always selected, because the observations are very time consuming. 

• The total number of regions in dark clouds, for which OH Zeeman observations are published, is still 
small. 

The best study of OH Zeeman splitting in dark clouds we are aware of is by Crutcher et al. (1993). 
They selected 4 cores in the Taurus dark cloud complex, 2 in the Libra complex, 2 in p Oph, 1 in the Orion 
molecular ridge, 1 position in L889, and the core of B1 (Barnard 1), in the Perseus region. The only certain 
detection is in the cloud Bl. For the other regions, the weighted average value of the field is +2.7 ± 1.5 /iG in 
Taurus, —2.1 ± 2.8pG in Libra, +6.8 ± 2.5pG' in pOph, —0.6 ± 2.1 pG in L889, and —4.7 ± 3.5pG in L1647. 

Two well known regions of very intense magnetic field are Orion A and Orion B. OH Zeeman splitting 
in Orion A revealed a magnetic field strength of B = —125pG (Troland, Crutcher, & Kazes 1986), and 
B = +38pG toward Orion B (Crutcher, & Kazes 1983). 

Crutcher et al. (1996) have reported the first attempt to measure CN Zeeman splitting. They observed 
towards the cores OMC-N4 and S106-CN. They did not detect the magnetic field, and their upper limits 
are several times smaller than the value expected if the observed motions were sub-Alfvenic. 

An average field of +9pG was found in Cas A by Heiles and Stevens (1986). 

The OH Zeeman splitting should probe regions with n = 10 3 — 10 4 cm” 3 , while the CN Zeeman 
splitting, should probe regions of n « 10 6 cm” 3 . 

Note that all works are biased toward regions of strong field, because only regions considered to be 
favorable for magnetic field detections are selected for Zeeman splitting observations. Despite that, in many 
cases the field is not detected, and upper limits are quite stringent. 

In our model A we find that, if the mean field strength is ( B 2 ) 1 / 2 = 5pG, and the mean density 
(n) = 10 3 cm” 3 , at the time t = 2.0 td yn , then 20% of the total mass is in dense clumps, with density ten 
times larger than the mean, n « 10 4 cm” 3 , and field strength five times, or more, larger than the mean, 

B > 25pG. Therefore, even if field strengths of about 30pG are detected sometimes in dense cores, the 
mean Alfven velocity in the molecular cloud may be just comparable to the sound speed, (A4) 1 / 2 ~ Cs- 

On the basis of the observational results by Crutcher et al. (1993, 1996), it may be concluded that 
model A is in consistent with the observational estimates of magnetic field strength in dark clouds. The 
particular values of the magnetic field used in the comparison here should not be taken too literally; the 
small scale field strengths could be larger than estimated by the Zeeman effect if the magnetic field is 
tangled. 

It is difficult to envisage how the measurements could be consistent with model B, however, since 
the magnetic field in a sub-Alfvenic model by definition has an energy that exceeds the kinetic energy. 
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As demonstrated by Galsgaard & Nordlund (1996), a magnetically dominated plasma is able to quickly 
dissipate structural complexity, independent of the value of the resistivity. Thus, a sub-Alfvenic field 
could not remain strongly tangled, and hence could not avoid detection. Using the same argument, we 
expect those cores where a strong (sub-Alfvenic) field has indeed been observed to have a relatively simple 
magnetic field structure. 

Model A can be used for the description of a typical molecular cloud with linear size of about 5 pc and 
(n) = 500 cm -3 . At the time t = 1.0 td yn , the relation B — n is then: 


B = 5 nG{- 




0.7 


(9) 


' 500cm 

Since the exponent is > 0.5, most of the dense cores are found to have magnetic pressure larger than 
thermal pressure, at early times. This is illustrated in Fig. 8, where we have plotted the Alfven velocity 
versus the gas density, after one dynamical time, for experiments A and B. In model A, one can see that, 
on average, regions with density larger than the mean are characterized by an Alfven velocity larger than 
the sound velocity, and therefore by a magnetic pressure larger than the thermal pressure. 


The lowering of the exponent in the B — n power law with time means that later on, the dominance 
of the magnetic field in these cores tends to be reduced. Although the flow is random and approximately 
isotropic, the kinematic alignment of B with v makes dense cores accrete mass along B at an increased rate. 
Therefore, the accretion of mass around dense cores, embedded in a random flow with (3 « 1, is such that, 
while magnetic pressure becomes dominant over thermal pressure during the initial phase of fragmentation, 
later on magnetic pressure decreases to the level of thermal pressure, even in the absence of gravity, and on 
a timescale that is competitive with ambipolar diffusion. This mechanism could be relevant for the process 
of star formation. 


The main point we want to clarify in this paper is that the fact of finding some cloud cores, with 
sub-Alfvenic velocity dispersions and with magnetic pressure larger than thermal pressure, does not 
necessarily mean that the dynamics of molecular clouds is dominated by MHD waves; those cores may be 
formed, in a few million years, in a supersonic and super-Alfvenic flow, only marginally affected by the 
magnetic field (model A). The dynamics becomes strongly affected by the magnetic field only in some very 
dense regions, on small scales, and at early times. 

We stress that a key theoretical ingredient to the interpretation of OH Zeeman splitting data is the 
B — n relation. The fact that the exponent of the relation is > 0.5 for more than two dynamical times is 
the reason why sub-Alfvenic cores can be found in the experiment. 


5. Conclusions 


In this work we have shown that: 

• Supersonic motions in model A are relatively long-lived, with decay times comparable to the estimated 
life-time of molecular clouds. Even in the absence of an energy source, supersonic and super-Alfvenic 
motions can persist as such, on a scale of 20 pc, for about 10 7 years. 

• Supersonic sub-Alfvenic motions (model B) dissipate their energy almost as fast as super-Alfvenic 
motions. The original motivation for the MHD wave model is therefore absent. 
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• Random supersonic motions produce a very intermittent probability distribution of the magnetic 
energy, with an exponential tail. Since the distribution is so intermittent, one could easily detect a 
held strength that is several times in excess of the mean strength. 

• Dense cores, with magnetic pressure larger than thermal pressure, and velocity dispersions smaller 
than va , are found as the result of the evolution of supersonic and super-Alfvenic hows (model A). 

• A power law statistical B — n relation is generated by supersonic motions, in model A, but not by 
MHD waves, in model B. The exponent of the relation is > 0.5 for about two dynamical times, which 
allows for the existence of cores with va > Cg. 

• The exponent in the B — n relation decreases with time, because magnetic held lines are stretched and 
partially aligned with the how. The statistical importance of the magnetic pressure in dense cores is 
thus expected to decrease with time, even in the absence of gravity and ambipolar diffusion. 

• The topology of the density held generated by model A is mainly filamentary and clumpy. This is 
reminiscent of the morphology of molecular clouds. On the other hand, the density held produced by 
MHD waves in model B is mainly structured in sheets, perpendicular to the magnetic held. 

• The statistical properties of the density held generated by random supersonic motions (as in model 
A) are in agreement with the statistical properties of the gas density in dark clouds, as inferred from 
stellar extinction determinations. Model B does not give such a good agreement. 

The comparison between model A and model B shows that model A (super-Alfvenic motions) provides 
a nice understanding of the dynamics and structure of molecular clouds. Moreover, it is consistent with the 
OH Zeeman detections. 

We conclude that the supersonic motions observed in dark clouds, on a scale of a few parsecs to several 
parsecs, can be understood, without any major problem, as super-Alfvenic motions, although sub-Alfvenic 
dense cores certainly exist. 
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Fig. 3.— Distribution of magnetic energy in units of the mean gas pressure, after one dynamical time. The 
dotted vertical line marks the initial mean value of the magnetic energy, and the dashed line the mean value 
after one dynamical time. 
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Fig. 4.— The relation B — n, after one dynamical time. The magnetic field is expressed in the numerical 
units, that is as an Alfven velocity (defined with the mean gas density), in units of the sound speed. The 
1 — a ‘error’ bars show the dispersion of values of B around the mean for each bin. The size of the bins is 
chosen in order to have the same number of measurements inside each bin. 
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Fig. 5.— 3-D representation of the time evolution of the density field. Time evolves from left to right. The 
upper row of panels show model A, while the lower row shows model B. 
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Fig. 6.— Histograms of the cosine of the angle between v and B. In experiment A there is a partial 
alignement, while in experiment B the two fields are mainly perpendicular to each other. 
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Fig. 7. Relation between the dispersion of the extinction and the mean extinction, on a regular grid 
superposed of the projection of the density field. The value of the extinction are measured at random 
positions, simulating the random position of a star behind the cloud. 
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Fig. 8.— The Alfven velocity, in units of the sound velocity, versus the gas density. The dashed line is the 
B — n relation, that is to say the Alfven velocity defined with the mean density, instead of with the local 
density. In Experiment A, the magnetic pressure is on average larger than the gas pressure, in regions with 
gas density larger than the mean. 





